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ABSTRACT 

We report on a multi-wavelength study of the recently discovered X-ray transient 
XMMU J004243.6H-412519 in M31, based on data collected with Swift and the 1.8-m Coper- 
nico Telescope at Cima Ekar in Asiago (Italy) between 2012 February and August. Unde- 
tected in all previous observations, in 2012 January XMMU J004243.6H-412519 suddenly 
turned on, showing powerful X-ray emission with a luminosity of ^10'^^ erg s~^ (assum- 
ing a distance of 780 kpc). In the following weeks, it reached a luminosity higher than '^10'^^ 
erg s~^, in the typical range of ultraluminous X-ray sources (ULXs). For at least '^40 days 
the source luminosity remained fairly constant, then it faded below «10'^^ erg s^^ in the fol- 
lowing ~^200 days. The source spectrum, which can be well described by multi-color disk 
blackbody model, progressively softened during the decay (the temperature changed from 
kT - 0.9 keV to -0.4 keV). No emission from XMMU J004243.6H-412519 was detected 
down to 22 mag in the optical band and to 23-24 mag in the near ultraviolet. We compare the 
properties of XMMU J004243 .6h-4 1 25 1 9 with those of other known ULXs and Galactic black 
hole transients, finding more similarities with the latter. 
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1 INTRODUCTION 



Ultra luminous X-ray sources (UL Xs; iLong & van Spevbroeckl 
ll983l : lHelfanjl984lFab"biandll989l) are observed in nearby galax- 
ies as point-like, off-nuclear X-ray sources with isotropic lumi- 
nosity larger than the Eddington limit for spherical accretion of 
fully ionised hydrogen onto a ~1O-M0 compact object (~10^' 
erg s~^). The nature of these objects remains a major astrophys- 
ical puzzle. Some ULXs have been identified with X-ray bright su- 
pernovae, other with background active galactic nuclei, but most 
of them are s uspected to be accre ting black holes (BHs) of stel- 
lar origin (see iFeng & Sorij|20T l] for a recent review). However, 
their luminosities are typically too high to be explained by accre- 
tion processes on stellar-mass black holes without violating stan- 
dard Eddington limit arguments. Three main different possibilities 
have been proposed to explain the nature of ULXs. (1) They are 
intermediate-mass BHs (with mass in excess of 100 Mq) in accre- 
tion ( Colbert & Mushotzky 1999|); a strong case for the existence of 
such objects is HL X- 1 in ESQ 243-49, for which the lumi nosity ex- 
ceeds 10"^ erg s"^ jFarrell et alj2009l : lDavis et alj201 ih . (2) They 
are accreting stellar-mass BHs (~5-20 Mq) with substantial beam- 
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ing of their X-ray e mission I King et a l]l200lh or emitting above 
their Eddington limit (Begelma nl l2002l ;l iBegelman. King & Pringld 
[2006) . (3) They are massive BHs (~30-80 Mq) formed in low- 
metallicity environments from massive progenitors and accreting at 
or sli g htly above their Eddington limit /Mapelli. Colpi & Zampieril 



|2009l ; IZampieri & Ro berts 2009; Belczynski e t al.| l2010h . Othei- 



possibilities are that the anomalously high luminosity of these ob- 
jects is due to a combination of these factors, or that the population 
of ULXs is composed by different kind of sources. 

Observations of ULXs' counterparts in energy bands other 
than X-ray have proven difficult because of their extragalac- 
tic nature. In X-rays, ULXs have been extensively observed 
during the last decade with high-spectral-resolution satellites 
such as Chandra and XMM-Newton. High-quality spectra taken 
with such satellites have shown that many bright persistent 
sources emit a spectru m commonly well fit by a Comptoniza- 
tion model plus a disc i Stobbart. Roberts & Wilmsll2006l ; [Roberts! 
|2007| ; iGladstone. Roberts & Pond |2009|) . These spectral proper- 
ties appear consistent with super-Eddington or marginally super- 
Eddington accretion onto a stellar-mass o r a more massive 
BH (e.g. [Gladstone. Roberts & Doiid [20091 ; [Pintore & Zampieril 
[2OI2I) . Quasi-periodic oscillations (QPOs) similar to the ones 
observed in the emission of Galactic BH binaries (see e.g. 
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ICasella. Belloni & Stellall2005[) h ave been detected in X-rays in 
a few ULXs jMucciarelli et alj |2006|. IStrohmaver et^ |2007| ; 



IStrohmaver & Mushotzkvll200ll2009l : lFeng. Rao & Kaarefcoiof) . 
Some efforts have been made to mak e use of scaling arguments 
valid for Galactic BH binaries (see e.g. lShaposhnikov & TitarchukI 
l2009h in order to constrain the masses of ULXs from characteristic 
time-scales. Results are, however, highly uncertain, as all these esti- 
mates are based upon tentative identifications of the timing features 
and the use of scaling l aws that are known to hold only for a limited 
number of objects (e.g. lMiddleton. Sutton & Robertsll201lh . 

Despite the substantial body of work that has been done 
in the last decade, transient ULXs in particular still repre- 
sent a poorly know population. Only few ULXs with tran- 
sient behaviour have been detected so far and only a hand- 
ful of them have been studied in detail. Examples of well- 
studied genuine ULX transients are CX0M31 J004253. 1+4 11422 
in M31 (hereafter M31 ULX-1), CXOU 133705. 1-295207 inM83 
and CXOUJ132 5 18.2^30304 in NG C5128 l lKaur et al.l [201I : 
ISoria et all [20121 ; Isivakoff et alj I2OO8I) . In 2012 January a new 
X-ray source with a luminosity of ^10'^* erg s"^ was dis- 
covered b y XMM -Newt on in M 31 (XMMU J004243.6+412519; 
iHenze. Pi etsch & Haberl. ii2012bl) . Seven days after its discovery, it 
reached a luminosity of ~ 2 x 10 ^^ erg s~^, which made it the sec - 
ond most luminous ULX in M31 faenze. Pietsch & Haberll2012ch . 
During the outburst the source was repeatedly observed by Swift 
and other instruments. Here we present the analysis of the sev- 
eral Swift XRT and UVOT observations performed between 2012 
March and August with the aim to explore the nature of this source. 
We also present the results of deep optical observations taken on 
2012 July 18 with the 1.8-m Copernico Telescope at Cima Ekar 
(Asiago, Italy) and the analysis of serendipitous observations of 
the source region collected in recent years by Swift and Chandra. 

In Section [2] we describe the Swift observations used in our 
study and we present the results of the spectral and timing analy- 
sis of the X-ray data in Section [3] In Section |4] we report on the 
SwiftfUVOT and Copernico optical and ultraviolet observations of 
the source. In Section [5] we present the upper limits on the pre- 
outburst X-ray flux of XMMU J004243.6+4125I9 obtained from 
our inspection of Chandra and Swift archival observations. Discus- 
sion follows in Section[6l 



2 Swift OBSERVATIONS AND DATA REDUCTION 

The Swift payload includes a wide-field instru ment, the coded- 
mask gamma-ray Burst Alert Telescope (BAT; iBarthelmv et al.l 
I2OO5I) , and two narrow-fi eld instruments, the X-Ray Telescope 
(XRT; lBurrows et al.ll2005l) and the Ultra- Violet/Optical Telescope 
(UVOT; lRoming et alj20^ . In this work we made use of only the 
narrow-field instruments data. 

The XRT uses a front-illuminated CCD detector sensitive to 
photons between 0.2 and 10 keV. Two main readout modes are 
available; photon counting (PC) and windowed timing (WT). PC 
mode provides two dimensional imaging information and a 2.507-s 
time resolution; in WT mode only one-dimensional imaging is pre- 
served, achieving a time resolution of 1.766 ms. The UVOT is a 
30-cm modified Ritchey-Chretien reflector using a microchannel- 
intensified CCD detector which operates in photon counting mode. 
A filter wheel accommodates a set of optical and ultraviolet (UV) 
filters and the wavelength range is 1700-6000 A. The data were 



processed and filtered with standard procedures and quality cutfl 
using FTOOLS tasks in the HEASOFT software package (v. 6. 12) and 
the calibration files in the 2012-02-06 CALDB release. 

Foll owing the discovery of XMMU J0 04243.6+4I25 19 (2012 
January; iHenze. Pietsch & Haberll l20I2lj) . the source was ob- 
served by Swift ten times in about three weeks, until it came 
out of visibility in 2012 March. At the end of 2012 May, 
XMMU J004243.6+4125I9 became visible again for Swift and the 
monitoring w as resumed with ma ny further pointings, mostly off- 
axis (see also iHenze et aljr2012l) . Around mid 2012 August, the 
source flux became too low for the typical sensitivity of a ~2-ks 
XRT snapshot. For this reason, after the visibility gap we consider 
only the sixteen observations taken up to 2012 September 01. A 
summary of the observations used in this work is given in Table[T] 



3 X-RAY DATA ANALYSIS 

We extracted the PC source events from a circle with a radius of 
20 pixels (1 XRT pixel corresponds to about 2" 36) and the WT 
data from a 40 x 40 pixels box along the image strip. To esti- 
mate the background, we extracted PC and WT events from re- 
gions far from the position of XMMU J004243.6-I-4125I9. The an- 
cillary response files (arf) were generated with XRTMKARF, and 
they account for different extraction regions, vignetting and point 
spread function corrections. We used the latest available spectral 
redistribution matrix (rmf) in CALDB. The spectral channels were 
grouped so as to have bins with a minimum number of 20 pho- 
tons. The spectral analysis was performed with the XSPEC 12.7 
fit ting package; (A rnaud 1996); the abundances used are those 
of lAnders & Grevessel lll989|) and photoelectric absorption cross- 
sections are from Balucinska-Church & McCammonI ( Il992l) . 

Initially we focus our analysis on the 2012 February-March 
data (see Table [T). For a preliminary look at the data, we fit all 
spectra simultaneously (in the 0.5-10 keV energy range) with the 
hydrogen column densit y tied between all observations using a 



hydrogen column densit y tied between all observations using a 
simple power law (e.g. iHenze. Pietsch & Haberll l2012d) . While 
this simultaneous modelling yields a rather high reduced of 
1.37 for 553 degrees of freedom (dof), the test shows that, as 
can be seen in Fig. [T] the spectra of XMMU J004243. 6-1-4 125 19 
are all very similar in this se t of XRT observations (see also 
iHenze. Pietsch & Haberill2012iJ) . 

Thus, in order to achieve better statistics and higher signal- 
to-noise ratio, we merged the data from all the observations and 
accumulated combined PC and WT spectra. In the following, we 
concentrate on the PC total spectrum (24.9 ks, about 7900 counts, 
99.4% of which are attributable to XMMU J004243. 6-^4 125 19), 
because of the intrinsically higher signal-to-noise ratio of the PC 
data with respect to the WT ones. The combined spectrum was fit 
adopting several phe nomenological models frequently used in liter- 
ature for ULXs (e.g. lRobertsl2007l;lFeng & Sorial201lh : a DISKBB, 
a POWERLAW, a DISKBB+POWERLAW, and a DISKBB+COMPTT, 
all corrected for interstellar absorption. All these models but the 
power law provide statistically good fits. However, the improve- 
ment given by the combination of two models is not signifi- 
cant in comparison with a single DISKBB model (x^ = 1-07 
for 231 dof against, for instance, = 1-03 for 229 dof for 
the DISKBB-l-POWERLAW). We also tried to fit more sophisti- 
cated Comptonization models (COMPTT and SIMPL in XSPEC) 



^ See' http://swift.gsfc.nasa.gov/docs/swift/analysis/| for more details. 
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Table 1. Swift observations used for this work. 



Observation XRT UVOT Start / end time (UT) Exposure 

mode filter" (YYYY-MM-DD hh-mm-ss) (ks) 
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" u: central wavelength 3465 A, FWHM 785 A; uvwl: central wavelength 2600 A, FWHM 
693 A; uvm2: central wavelength 2246 A, FWHM 498 A; uvw2: central wavelength 1928 A, 
FWHM 657 A. 



Table 2. Spectral analysis of XMMU J004243.6-^412519 (2012 February-March, PC data). Errors are at a Icr confidence level for a single parameter of interest. 



Model" Wh Ti -Ecut/break kTin r2 / Norm.'' Observed fliix'^ Luminosity'^ xi (dof) 

(10^1 cm-2) (keV) (keV) (10"" erg cm-^ s^i) (lO^^ergs-1) 



PHABS*POWERLAW 6.7 ± 0.2 2.91 ± 0.04 - - - 1.12 ±0.02 2.33 ±0.09 1.46(231) 

PHABS*DISKBB 2.9 ±0.1 - - 0.86 ± 0.02 1.4 ±0.1 1.04 ±0.02 l-O'ito^^l 1.07(231) 

PHABS*cuTOFFPL 3.5 ± 0.3 0.7^°^ ^-^-oi - - ^-OSto.oi l-15io:o6 1.03 (230) 

PHABS*BKNPOWER 4.5 ± 0.3 ^-itgl ^-^toA - -^-^-o.! 1.10 ±0.02 l-44iS;o8 1.00(229) 

PHABS*(POWERLAW-HDISKBB) 3.2^02 l-^t^s - ^■'^^-0 02 1.8 ± 0.2 1.10 ±0.02 1.121°°^ 1.03 (229) 

PHABS*(CUTOFFPL-fDISKBB) 9^^ 1.4 ± 0.5 l-7to,l 0.12t_°g°l l.lt^J X 10^ l-07to.ol l-15±S;o6 1.00 (228) 



" XSPEC model. 

^ DISKBB normalisation; (7?.[krn]/£>[10 kpc])-^ cos 6,9 — coiTesponding to a face-on disk. 

^ In the 0.5-10 keV energy range; for the luminosity we assumed a distance to M3 1 of 780 kpc lHolland[l998l ; IStanek & Garnavichlll998l) . 



to the data, but the count statistics of the spectrum is not suffi- 
cient to provide good-enough constraints on the model parame- 
ters. Hence, it seems appropriate to describe the pre-gap PC spec- 
trum of XMMUJ004243. 6-1-412519 in terms of a simple DISKBB 
model. The best fitting parameters show that the disc is rela- 
tively hot (0.86 ± 0.02 keV) and the luminosit y is quite high 
(~1.1 X lO^'^ for a distance of 780 kpc; lHolland|[T99a ; 

IStanek & Garnavicl]|l998h . The resuUs are summarised in Table[2] 



Similarly, for the second batch of data (collected starting 
from 2012 May; see Table [TJ we extracted a cumulative PC 
spectrum (23.7 ks, about 1600 counts, 99.2% of which are at- 
tributable to XMMUJ004243.6-I-412519). We tested the same 



single-componen(3 spectral models reported in Table [2] obtaining 
the parameters shown in Table [3] In general, they appear to have 
values of the column density consistent with those of the pre-gap 
state and significantly lower fluxes/luminosities. In particular the 
post-gap average spectrum is well described by a DISKBB com- 
ponent with an (average) temperature (~0.6 keV) smaller than 
that of the pre-gap combined spectrum. In Fig. |2] we show the 
PHABS*DISKBB model fit to the pre- and post-visibility gap com- 
bined PC spectra. The softening suggested by the spectral parame- 
ters and Fig. [2] is apparent from the contour plots shown in Fig. [3] 



^ Owing to the lower count statistics, this time we did not consider more 
complicated (two-component) models. 
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Table 3. Spectral analysis of XMMU J004243.6+412519 (2012 May-August, PC data). Errors are at a Icr confidence level for a single parameter of interest. 



Model" JVh Ti Bcut kTi„ r2/Norm.'' Observed flux"^ Luminosity"^ X? (dof) 

(10^1 cm-2) (keV) (keV) (lO^i^ erg cm^^ j-l) (lO^^ergs-1) 

PHABS*POWERLAW 6.9 ± 0.5 3.6 ±0.1 - - - 2.04 ±0.07 0-7lioJs 1.45 (66) 

PHABS*DISKBB 2.7 ± 0.3 - - 0.62 ± 0.02 1.2 ±0.2 1.95 ± 0.06 0.22 ± 0.01 0.98 (66) 

PHABS*cuTOFFPL l-Slas ^l-Oiae 0.6 ±0.1 - - 1.95 ±0.06 0-17iao2 0.98 (65) 

PHABS*BKNPOWER 2.2tJ5;^ '^'^'^0^3 - ^.e+jj;^ 2.0lt"a°^j 0.29 ±0.05 0.99(64) 

" XSPEC model. 

^ DISKBB normalisation: (i?[km]/i?[10 kpc])-^ cos &,8 — corresponding to a face-on disk. 

^ In the 0.5-10 keV energy range; for the luminosity we assumed a distance to M3 1 of 780 kpc lHolland[l998l : IStanek & Garnavichlll998l) . 
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Figure 1. Simultaneous modelling of all the 2012 February-March PC and 
WT spectra using an absorbed power-law model. Bottom panel: the residu- 
als of the fit (in units of standard deviations). (See the electronic journal for 
a colour version of this figure.) 




nH (10=^ cm-') 

Figure 3. Contour plots of temperature versus column density (adopting the 
PHABS*DISKBB model) for the pre-gap (black/red/green upper contours) 
and post-gap (blue/cyan/violet lower contours) PC data. The crosses indi- 
cate the best fits. Lines mark 1, 2 and 3(T confidence levels. 
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Figure 2. Comparison of the pre- (black circles) and post-visibility gap 
(blue triangles) cumulative PC spectra for the PHABS*DISKBB model (see 
Section |3] for details). Bottom panel: the residuals of the fit (in units of 
standard deviations). 

In order to study the source flux evolution over time, we 
fit afl the PC spectra from the individual observations or from 
small groups of observations. When the source flux dropped below 
2 X 10^^^ erg cm~^ s~^, around mid 2012 July, in order to ac- 
cumulate sufficient statistics for meaningful spectral fits, we com- 
bined data from a few contiguous observations; namely, we ob- 
tained a spectrum from observations 00035336058-61 and another 
one from observations 00035336062-64. We fit all the spectra si- 



multaneously adopting the PHABS*DISKBB model with the hydro- 
gen column density (which is consistent with a single value in both 
the pre- and post-gap cumulative spectra) tied between observa- 
tions (xS = 1.10 for 676 dof). We plot the resulting long-term light 
curve and the inferred temperatures kT in Fig.|4] We show both the 
absorbed and unabsorbed fluxes. As can be seen from Fig. |4] the 
effect of the interstellar absorption is larger on the softer spectra 
from the post-gap observations. Progressive spectral softening and 
flux decay are evident during the post-gap observations. The flux 
decreased by a factor of ~5 over ^^70 days (by a factor of ~10 
over ~150 days with respect to the pre-gap flux), while the disk 
temperature changed from kT ^ Q.7 keV to ~0.4 keV. 

Although, owing to the long visibility gap, the available data 
do not allow us to perform an accurate modelling of the decay 
shape, we tried a number of simple models to fit the light curve 
of XMMUJ004243.6-I-412519. We fixed t = Q si the time the 
source was observed for the first time to exceed the U LX thresh- 
old (MJD 55947.5 h iHenze. Pietsch & Haberll [20123) . An expo- 
nential function of the form F{t) = Aexp(— f/r) gives a rather 
poor fit (x? ~ 11.58 for 10 dof). The best-fitting parameters for 
the observed (absorbed) flux are A = (2.00 ± 0.05) x 10"" 
erg cm"^ s"^ and e-folding time r — (63.9 ± 1.3) d. 

Among the other models tested, a broken-power-law model 
provides the better fits for the evolution of both the absorbed 
(X^ = 4.33 for 8 dof) and unabsorbed fluxes (x2 = 3.65 for 8 
dof). Assuming as t = MJD 55947.51, for the observed flux, 
the break occurs at (111.4 ± 6.4) d, when the index changes from 
ai = -0.41 ± 0.18 to Qi = -4.03 ± 0.19; the flux at the 
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Figure 4. The top panel sliows time evolutions of ttie absorbed (black cii- 
cles) and unabsorbed (red triangles) fluxes in the 0.5-10 keV energy range 
for the PHABS*DISKBB model (see Section[3]for details); the down-arrows 
indicate upper limits at the 3cr confidence level. The broken-power-law 
models describing the decays are also plotted (black solid line for the ob- 
served flux and red dashed line for the unabsorbed flux). Bottom panel: evo- 
lution of the characteristic temperature of the DISKBB model inferred from 
the spectral fitting. We assumed as t = the dat e the source was observed 
for the first time to exceed the ULX threshold iHenze. Pietsch & Haberll 
l2012cl) . 



To improve the statistics, we produced a single P DS averaging 
all the WT data. The PDS was normalised according to lLeahv et al.l 
( Il983h , so that powers due to Poissonian counting noise have an 
average value of 2. Since PDS created from data taken in WT mode 
show a drop-off at high frequencie^ and since also the averaged 
PDS from all the WT observations did not show any significant 
feature, we rebinned our data at a Nyquist frequency of ~17 Hz 
(rebin factor of 16) using data stretches 925-s long to obtain a new 
average PDS. The PDS is well fit by a constant component equal 
to 2.286 ± 0.006 (reasonably describing the Poissonian noise) with 
a best-fitting xi^ = 169.56 for 152 dof. Hence we conclude that no 
significant variability is detected from the source. It is not possible, 
however, to exclude that the emission from the source has a certain 
level of variability hidden by photon counting statistics or on time- 
scales not accessible to our data (<10"'' Hz). 

Apart from the flux decay on the scale of weeks, the emission 
from XMMUJ004243. 6-1-412519 does not show signs of strong 
aperiodic variability either. There is only some evidence of mod- 
erate variability on the 1-ks-scale, with a rms variance of (20 ± 5) 
percent in the WT data. This is probably related to some level of 
variability of the local astrophysical background and could account 
for the fact that the Poissonian noise level observed in the PDS is 
slightly higher than the expected value (2 with the Leahy normal- 
isation). For the pre-gap PC light curves (a similar analysis of the 
post-gap data is hampered by the low count rates) the 3a upper lim- 
its on the rms variability range from approximately 12 percent to 
20 percent. 



break time is (6.9 ± 1.5) x 10 erg cm ^ s ^. For the unab- 
sorbed flux, the best-fitting parameters are: Qo,i — —0.26 ± 0.15, 
ao,2 = -3.06 ± 0.19 and (104.5 ± 6.4) d for the break epoch 
(with an unabsorbed flux of (1.1 ± 0.2) x 10~^^ erg cm~^ s~^). 
The steeper decay of the absorbed flux is due to the fact that the ab- 
sorption affects the softer (later) spectra more. We stress that the fit 
parameters depend on the assumed time origin. The uncertainties 
introduced by making different (reasonable) assumptions can be 
larger than the statistical errors reported here. For instance, assum- 
ing t = at MJD 55942 (the day after t he discovery of the source 
at a lu minosity of ~ 2 x 10^* erg s~^: [Henze. Pietsch & Haberll 
l2012bl) . we find for the observed flux: a[ = -0.19 ± 0.12, 
Qa = -4.18 ± 0.19, (110.3 ± 3.6) d for the break epoch and 
flux at the break time of (8.8 ± 1.0) x 10"^^ erg cm"^ s~\ 

In the last observations (segments 00035336065, 
00035336066 and 00035336067), XMMU 1004243.6-^412519 
was not detected. The 3a upper limits on the XRT count rate 
derived from the deeper observations (6065 and 6067) are of 
~0.1 counts s^^ (0.3-10 keV, following lKraft. Burrows & Nousekl 
Il99 ih . Assuming the spectrum of the closest observations (6062- 
6064, kT ~ 0.4 keV), this corresponds to upper limits on the 
observed flux of 

erg s^^ on the luminosity (for 780 kpc) 

For the timing analysis we concentrate on the data in WT 
mode (total exposure: 16.0 ks; see Table [T), since their time res- 
olution of -^1.7 ms (corresponding to a Nyquist frequency of 
^280 Hz) makes it possible to search for fast time variability. 
For each observation we computed a power density spectrum 
(PDS) in the energy band 0.3-10 keV by using intervals up to 
~925-s long and averaging the individual spectra for each obser- 
vation. The PDS did not reveal any significant variability from 
XMMU J004243.6-I-4125I9: neither broad-band components nor 
narrow features were detected; no significant evolution of the PDS 
in time could be detected as well. 



-4 X 10"^^ ergcm"^ s"^ and of ~6 x 10^^ 



4 OPTICAL AND ULTRAVIOLET OBSERVATIONS 

The Swift/UVOT observed XMMU J004243.6-F412519 simultane- 
ously with the XRT. The data were taken with the u, uvwl, uvm2 
and uvw2 filters (see Table [T]l. The analysis was performed on 
the individual and stacked (for each filter) images with the UVOT- 
SOURCE task, which calculates the magnitude through aperture 
photometry within a circular region (we used a 3-arcsec radius) 
and applies specific corrections due to the detector characteristics. 

No source was detected at the position of 
XMMUJ004243.6-(-412519 in any of the UVOT observations 
and filters, before or after the visibility gap. The 3a limits before 
the visibility gap are in the stacked images u > 23.2 mag (total 
exposure: 13.4 ks), uvm2 > 23.7 mag (total exposure: 5.6 
ks) and uvw2 > 24.5 (total exposure: 18.2 ks). After the gap, 
XMMUJ004243.6-I-4125I9 was observed only with the uvwl 
filter; the 3a limit from the stacked image is uvwl > 24.4 
(total exposure: 23 .2 ks ) All magn i tudes are in the AB system 
( IOke&Gunrilll983l) ; see lPoole et aj] (l2008h for more details on the 
UVOT photometric system and lBreeveld et al. l|) for the most 
updated zero-points and count rate to flux conversion factors. 

The above magnitudes have not been corrected for extinction. 
At the position of XMMU J004243.6-I-4I2519, the total line-of- 
sight optical extinction esti mated from background infrared emis - 
sion is Av ~ 0.05 mag dSchlegel. Finkbeiner & DavisI Il998h . 
while the X-ray fits, adopting the rel ation A'^h = 1-79 x lO^Mv 



cm ~ by IPredehl & Schmitd ( Il995h , yield higher values in the 
range Av ~ 1.7-5.0. Indicatively, Ay = 1 mag corresponds to 
Au ^ 1.8 mag, Anvwi — 2.2 mag , Auvm 2 — 3.2 mag, and 
Auvw2 ^ 2.8 mag jFitzpatrick & Mass j|2007h . 

^ This effect is mostly visible above ~ 50 Hz and is related to the read-out 
method of the XRT detector, see |http://www.swift.ac.uk/analysis/xrt| 
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Figure 5. V-band image of the field of XMMU J004243.6+412519 taken 
on 2012 July 18 at the 1.8-m Copernico Telescope at Cima Ekar in Asiago. 
North is up. East to the left. The green circle (7 "2, 3<7) is centred o n the 
Chandra position of XMMU J004243.6+412519 ifiamard et al.l2012l) . 



The field of XMMU J004243.6+412519 was also observed 
with the 1.8-m Copernico Telescope at Cima Ekar in Asiago (Italy) 
on 2012 July 18. Three images of 20 minutes were taken in both 
the V and B band filters. The data were reduced following standard 
prescriptions. After removal of the detector signature (bias and flat 
field corrections), the images were astrome trie ally calibrated per- 
forming a polynomial interpolation starting from the positions of 
the NOMAD star catalogue jZac harias et al. 200^ . The accuracy 
is 0'.'2. The three calibrated frames in each filter were then aver- 
aged and the resulting V band image is shown in Fig. [5] 

Instrumental magnitudes were measured on the images 
through the point-spread function (PSF) fitting technique. The pho- 
tometric calibration was perform ed using refere n ce sta rs from the 
catalogue of M31 compiled by iMagnier et alj jl992l) . homoge- 
neously distributed around the source position. The internal accu- 
racy of this calibration is 0. 1 mag in both bands. No source at the 
position of XMMUJ004243.6-I-412519 was detected in any of the 
two filters down to a limiting magnitude of 21.7 and 22.2 in the V 
and B band filters, respectively (see Fig.O. In fact, the background 
emission from M3 1 is highly variable inside the error box, so that 
the actual limiting magnitude varies from 21.5 to 21.9 in V and 
from 21.6 to 22.8 in B, depending on the position. 



5 PRE OUTBURST OBSERVATIONS 

As noted bv lHenze. Pietsch & Haberil ( I2012bl l3). no source compat- 
ible with the position of XMMU J004243.6-I-412519 was listed in 
any X-ray catalogue. We have searched the XMM-Newton, Chan- 
dra, and Swift public archives for possible previous bright states 
of XMMU J004243. 6-1-412519 in recent years, but the source was 
never detected. Since the XMM-Newton observations are rather 
sparse, in the following we summarise only the upper limits ob- 
tained from the much more intense coverage (from 1999 November 
to 2012 January) with Chandra and Swift. 
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Figure 6. Chandra 3 (7 upper limits on the 0.5-10 keV lumino sity (for a 
distance of 780 kpc: lHollanj|l998t Istanek & Gamavichlll998h . (See the 
electronic journal for a colour version of this figure.) 



were carried out with the HRC-I instrument I Murrav et alj|200(3h . 
39 with the ACIS (I or Si lGarmire e"tai]|2003ir Typical exposures 
are for the HRC-I ~ 1-5 ks in the 1999-201 1 observations and ~20 
ks in the more recent ones, and ~5 ks for the ACIS pointings; the 
deepest upper limit (see Fig.|6j was obtained from a 38-ks ACIS-I 
observation carried out on 2001 October 05 (obs. ID: 1575; MJD 
52187). 

For each observation, an upper limit on the count rate from 
XMMU J004243.6-I-412519 was computed using the CIAO tool 
APRATES and taking into account the point-spread function frac- 
tion in the apertures within which the counts were extracted. In 
order to convert limits on the count rates from the different de- 
tectors into upper limits on the unabsorbed flux, we used the 
NASA/HEASARC PIMMS tool assuming a soft power-law spec- 
trum with photon index F ~ 3 and absorption A*'h ~ 7 x 10^^ 
cm^^. The corresponding upper limits on the luminosity (assum- 
ing a distance of 780 kpc) range from ~6 x 10''^ to 5 x 10^^ erg s~^ 
and are shown in Fig. [6] 



5.2 Swift 

In the period between 2006 September 01 and 2012 January 01, 
Swift serendipitously imaged with the XRT (in PC mode) the posi- 
tion of XMMU J004243.6-I-412519 119 times for a total exposure 
of ~370 ks. The spread of the observations can be seen in Fig. |7] 
while the exposure time for each year is given in Table|4] 

We examined each observation, but the source was never de- 
tected, nor it was detected in the total and yearly total images. For 
each year and for the total data set, we computed 3a upper h mits on 
the count rate (following iKraft. Burrows & Nousel3ll99ll) . These 
upper limits (Table |4ll can be directly compared with the average 
(and fairly constant) PC rate observed in 2012 February-March, 
(0.370 ± 0.007) counts s^^. In particular, the last Swift observa- 
tion performed before the discovery of XMMU J004243. 6-1-4 125 19 
(obs. ID: 00035336051; 2012 January 01) yields an upper limit of 
2.5 X 10^'^ counts s^^, implying for the source a flux increase of 
150 or more during 2012 January /February. 



5.1 Chandra 

The distribution of the 93 Chandra observations (from 1999 
November 30 to 2011 August 25) covering the field of 
XMMU J004243.6-I-412519 can be seen in Fig.|6l 54 observations 



6 DISCUSSION AND CONCLUSIONS 

XMMU J004243. 6-1-4 125 19 is a transient X-ray source in M31 that, 
at its maximum, reached luminosity in the ULX range. Undetected 
in all previous X-ray observations, in February 2012 it suddenly 
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Figure 7. Epochs of the SwiftlXKI serendipitous observations of the position of XMMUJ004243. 6+412519. 



Table 4. Pre-outburst upper limits from Swi^/XRT serendipitous observa- 
tions of the position of XMMU J004243.6+412519. 



Year Observations Total exposure Count rate upper Umit" 



(ks) (counts s ^) 



2006 2 9.9 0.0010 

2007 19 58.3 0.00038 

2008 23 65.3 0.00023 

2009 28 86.2 0.00028 

2010 27 98.5 0.00031 

2011 19 47.9 0.00066 
2012 1 4.1 0.0025 

2006-2012 TTg 3693 0.00016 



° At a confidence level of 3o". 



started to show powerful X-ray emission. After reaching the max- 
imum, the source luminosity remained fairly constant at ^lO^" 
erg s~^ for at least ~40 days, then it faded below ~10^*^ erg s^^ 
in the following ^^200 days. The decay, accompanied by a spec- 
tral softening, can be described by a broken power-law with a 
break time of ~100 days. No broad components nor narrow fea- 
tures (in the form of quasi-periodic oscillations) were detected in 
the power density spectra up to ~280 Hz. We searched also for op- 
tical and near-UV emission from XMMU J004243. 6+412519, but 
no no source was detected at its position down to a sensitivity limit 
of ~22 in the V and B bands, and of I'i-l'X in the near UV. 

In the following, we discuss the properties of 
XMMUJ004243.6+412519 only in the context of an accret- 
ing black hole in M31. The large absorption column and the lack 
of an optical counterpart exclude a foreground object. In terms 
of a background object, the best candidates are an active galactic 
nucleus or a tidal disruption event. However, no active galactic 
nucleus has been observed to vary in X-rays by more than two 
orders of magnitude, while a tidal disruption event is ruled out by 
the flux evol ution, that does not follow the characteristic decay of 
such events jBurrows et al.ll20lll) . 

Our spectral analysis of the Swift PC data showed that the 
energy spectrum of XMMU J004243. 6+412519 can be fit using 
a disc component, and that the combination of a disk component 
and a power law that is commonly used to describe both the spec- 
tra of some ULX s and of Galactic BH binaries is not statistically 
needed (see also iBamard. Garcia & Murray! l2012h . However we 
cannot exclude the possibility that an additional component is not 
detected because of the low counting statistics. Along its outburst 
XMMU J004243. 6+412519 showed only moderate and slow spec- 
tral changes in contrast wit h what is usually seen in transie nt BH 
binaries in outburst (see e.g. lMotta. Belloni & Homaiill2009l) . 

The transient ULXs observed so far appear to have rather dif- 



ferent properties and XMMU J004243. 6+4125 19 is not an excep- 
tion in this picture. The source presents similarities with other tran- 
sients, but it cannot be definitely associated to any of them. 

M3 1 ULX- 1 was observed to reach a peak luminosity of 
5 X 10'"' erg s^^ jKaur et al.l2012h followed by a decrement in flux 
with an e-folding time-scale of ~40 days (Middleton et al. 201^. 
A blue optical counterpart to M31 ULX- 1 was detected during the 
outburst, but not in quiescence, pointing towards emission from an 
irradiated accretion disc. This suggests a low mass companion star 
transferring matter via Roche lobe overflow. The spectral properties 
of M31 ULX-1 are roughly described by a single disc component 
around a 10- Mq BH with a spin of a ~ 0.4|j On the other hand, 
contrarily to what we observed in XMMU J004243.6+412519, the 
best description of the spectra is obtained with the addition of 
a second component which represents an optically thick medium 
(namely a low temperature corona described by a COMPTT) close 
to the BH. It can be associated to the photosphere of a wind ejected 
by the inner regions of the disc, that is expected to set in at super- 
Eddington accretion rate. Indeed, the softer component of the spec- 
trum can be also described by a slim disc model, indicating that 
advection may be important. The disc component becomes more 
important when the total luminosity decreases, suggesting that the 
wind unveils progressively the inner regions of disc when the ac- 
cretion rate drops. 

Also M83 ULX-1 jSoria et al.ll2012h showed a long term evo- 
lution similar to XMMU J004243.6+4125I9 but, as in the case of 
M31 ULX-1, the energy spectrum is significantly different from 
XMMU J004243.6+412519. Undetected (Lx < lO^** erg s"^) be- 
fore 2011, M83 ULX-1 was discovered by Chandra and it was 
seen to reach a peak luminosity comparable to that of M3 1 ULX- 1 
(^^4 X 10"^^ erg s~^). Its spectral properties are modelled by a cold 
disc plus a power law without evidence of curvature at high en- 
ergy, opposite to what is commonly seen in ULXs. In addition, 
the source does not show hints of a decline of the emission for 
at least 100 days after discovery and no spectral transitions were 
observed. Again contrarily to XMMU J004243.6+412519, a blue 
optical counterpart is observed, but only during the outburst, point- 
ing to reprocessed emission from the outer accretion disc and the 
companion star, a red giant or AGB star with mass M < 4 Mq 
dSoria et al.ll2012h . 

Finally, the ULX in NGC5128, discovered by Chandra in the 
2009, was observed at a luminosity of (2-3) x 10^^ erg s^^. The 
outburst lasted for at least 70 days, but the source was no longer 
observed. It is an intriguing ULX because it experienced spectral 
state variations between a very high state, dominated by a power 



^ The dimensionless spin parameter a = cj / (GM^), where J and M are 
the angular momentum and mass of the black hole, respectively. 
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law with a photon index higher than 2.2, and a high/soft state, dom- 
inated by an accretion disc of '--^ 1 keV, consisten t with the behaviour 
of the Galactic BH XRBs jSivakoff et alfcOOSh . 

One possible interpretation of the spectra-timing properties of 
XMMUJ004243. 6+412519 can be given by making an hypothe- 
sis on the nature of the source. If it is a stellar-mass BH transient, 
we can assume that most of the outburst evolution of the source 
could not be observed due to its distance from the observer. In 
this scenario, the source would be visible from Earth only during 
its brightest phases, which for many BH transients (e.g. H 1743- 
322, XTEJ1650-7 54, GRO 11655^0) is encountered dti ring the 
soft spectral states jBelloni. Motta & Muiioz-Dariasll2oTll) . In the 
soft states the energy spectrum is strongly dominated by a soft disk 
component (that sometimes is the only component visible in the 
X-ray spectrum) while the fast time variability is usually consistent 
with zero. These properties are consistent with what we reported on 
XMMUJ004243.6+412519. The duration of this high-luminosity 
phase in a Galactic BH transient is variable depending on the source 
and on the properties of the single outburst. However, the aver- 
age length of such periods (few months) is absolutely consistent 
with the duration of the outburst of XMMU J004243.6+412519 
(see e.g. the case of the Galactic BH tra nsients XTE J 1550 -754, 
iKubota & Do"n3l2004l and GROI1655^0, lMotta et ani2012h . 

Assuming that at maximum luminosity 
XMMUJ004243.6+412519 was in a disc-dominated state 
and that it was radiating at a significant fraction of LecW (say 
~0.6), the mass of the BH would be ~12 Mq. For small inclina- 
tions (^45°), a similar value is obtained from the normalisation of 
the disc component obtained from the fit of the combined pre-gap 
PC spectrum, assuming that the inner disc radius is truncated at 
6 gravitational radii and that the dis c spectrum has a standard 
color co rrection factor (see eg. [2] in IZampieri & Roberts! [2OO9I ; 
see also iLorenzin & Zampierilhopgl and references therein). This 
is consistent with the hypothesis that XMMU J004243.6+412519 
could be indeed an accreting stellar-mass black-hole binary 
observed in its soft state. 

The upper limits in optical bands are sufficiently deep to 
place interesting constraints on the donor mass. Assuming no ex- 
tinction and a distance modulus of 24.47 mag for M31 (from 
the NASA/IPAC Extragalactic Database)Q the upper limit in the 
V band translates into an upper limit on the absolute magnitude 
Mv > —2.8. Taking binary evolution effects and X-ray irradia- 
tion into account, a stellar-mass BH accreting through Roche lobe 
overflow is consistent with this upper limit if t he donor is a main se- 
quenc e star of 8-10 A/0 or a giant of <8 A/0 i jPatruno & Zampieril 
|201(]|). In fact, for a donor below 5 Mq, the dis c is no longer stable 
I Dubus et al.|[l999l ; IPatruno & Zampierill2008h . in agreement with 
the transient nature of the source. 
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